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Highly  vibrationally  excited  w-propylbenzene  cations,  were  prepared  by  the  charge  transfer  reaction 

02+  +  C9H12  -*02  +  C9H121*  in  a  turbulent  ion  flow  tube.  The  subsequent  competition  between  fragmentation 
of  C9Hi2+*  into  C7H7+  +  C2H5  and  stabilization  in  collisions  with  N2  was  studied  at  temperatures  in  the 
range  423—603  K  and  at  pressures  between  15  and  200  Torr.  Most  of  the  C7H7+  is  the  aromatic  benzylium 
isomer,  while  the  fraction  of  the  minor  species,  seven-membered-ring  tropylium,  increases  with  T,  from  5  to 
20%.  Minor  fragments  are  CfTh^,  C7Hg+,  and  CgH9+.  Energy-transfer  step  sizes  (AE)  for  collisional  deactivation 
are  obtained  by  combining  the  stabilization  versus  fragmentation  ratios  measured  as  a  function  of  pressure  in 
this  study  with  fragmentation  rates  from  the  literature.  The  values  are  compared  with  related  information  for 
other  excited  molecular  ions  and  are  similar  to  those  of  their  neutral  analogues.  At  the  highest  temperatures, 
C9H12"*"  was  also  observed  to  pyrolyze  after  collisional  stabilization.  Employing  unimolecular  rate  theory,  the 
derived  rate  constants  for  thermal  dissociation  of  C9H12+  are  related  to  values  derived  from  the  specific  rate 
constants  k(E,J)  for  fragmentation.  Good  agreement  is  found  between  measured  and  predicted  pyrolysis  rate 
constants.  This  allows  us  to  confirm  the  dissociation  energy  of  into  C7H7+  (benzylium)  and  C2H5  as 

166.9  (±2.2)  kJ  mol-'  (at  0  K). 


1.  Introduction 

Charge  transfer  can  generate  chemically  activated  molecular 
ions  which  subsequently  undergo  a  variety  of  reactive  and 
nonreactive  processes.  Reactive  processes  include  dissociation 
at  high  excitation  energies  and  secondaiy  reactions.  The 
unreactive  or  energy-dissipative  processes  include  collisional 
stabilization  and  radiative  cooling  of  vibrationally  excited 
molecular  ions.  At  high  temperatures,  the  stabilized  molecules 
subsequently  undergo  thermal  dissociation. 

In  previous  work,  we  have  generated  highly  vibrationally 
excited  alkylbenzene  cations  by  charge  transfer  from  small  ions 
to  alkylbenzenes  in  a  high-temperature  flowing  afterglow 
(HTFA)  and  in  a  selected  ion  flow  tube  (SIFT).1-3  Branching 
ratios  were  determined  as  a  function  of  the  temperature,  and 
mechanisms  were  determined.  The  comparison  of  charge 
transfer  breakdown  curves  taken  at  about  1  Torr  to  data  from 
low-pressure  photodissociation  experiments  suggested  that  the 
buffer  gas  had  an  influence  by  either  collisionally  stabilizing 
and/or  thermally  dissociating  the  ions.  The  small  pressure 
variation  (at  most  0.2—1  Torr),  however,  did  not  allow  for  a 
systematic  analysis  of  pressure  effects.  The  subsequent  con¬ 
struction  of  a  high-pressure  turbulent  ion  flow  tube4’5  (TIFT), 
operating  at  buffer  gas  pressures  from  15  up  to  700  Torr,  allows 
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us  to  better  investigate  pressure  dependences.  Recent  improve¬ 
ments  to  the  TIFT  allow  measurements  over  a  broader  temper¬ 
ature  range  than  in  previous  studies.  Studies  of  this  type,  for 
example,  the  reaction  of  02+  with  ethylbenzene,  allowed  for 
an  investigation  of  the  competition  between  unimolecular 
fragmentation  of  ethylbenzene  cations  (forming  methyl  radicals 
and  C7H7+  isomers:  the  aromatic  benzylium  and  the  seven- 
membered-ring  tropylium  ions6)  and  collisional  stabilization  of 
the  excited  ethylbenzene  ions.7  Since  the  specific  rate  constants 
k(E)  for  fragmentation  are  known  from  separate  absolute  rate 
measurements,8,9  the  effective  rates  for  collisional  stabilization 
could  be  deduced.7  These  results  were  analyzed  in  a  simplified 
master  equation  approach10  leading  to  the  product  of  the  total 
collision  frequency,  Z,  and  average  step  size  for  deexcitation, 
(AE),  at  each  temperature.  (AE)  values  were  obtained  by 
identifying  Z  as  the  ion-induced  dipole  capture  rate  constant. 
The  (AE)  values  were  found  to  be  similar  to  those  for  the 
corresponding  neutral  alkylbenzenes. 

Little  information  exists  on  collisional  energy  transfer  from 
molecular  ions  which  are  highly  vibrationally  excited.  While 
relative  collision  efficiency  measurements  are  available,11-13  it 
is  difficult  to  put  them  on  an  absolute  scale  of  energy  transfer 
rates.  Since  our  02+  +  ethylbenzene  experiments7  have  led  to 
values  of  the  product  Z(AE),  it  appears  desirable  to  derive 
similar  information  also  for  other  excited  polyatomic  ions. 

This  is  the  aim  of  the  present  series  of  articles  in  which 
experiments  for  the  reactions  of  02+  with  n-propylbenzene  (this 
work)  and  with  n-butylbenzene14  are  described.  A  second  goal 
of  the  our  work  is  to  study  the  pyrolysis  of  the  stabilized 
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Stabilization  and  Dissociation  of  C9Hi2+  Ions 

alkylbenzene  ions  and  to  attempt  to  relate  specific  rate  constants 
k(E)  for  fragmentation  with  the  derived  thermal  dissociation  rate 
constants.  The  present  work,  involving  ^-propylbenzene  ions, 
presents  this  approach  for  the  first  time. 

The  chemistry  of  the  02+  +  «-propylbenzene  (C9H12)  system 
involves  numerous  processes:15 


02+  +  C9H12  — *  02  +  C9H12+* 

(1.1) 

C9H12+*  —  C7H7+  (Bz+)  +  C2H5  k2 

(1.2) 

C9H|2+*  —  C7H7+  (Tr+)  +  C2H5  Jfc3 

(1.3) 

C9H12+*  +  M  —  C9Hi2+  +  M  kA 

(1.4) 

C9H12+*-C9Hn++H  k5 

(1.5) 

c9h12+*  — *  C8H9+  +  CH3  k6 

(1.6) 

C9HI2+*-C6H6++C3H6  kn 

(1.7) 

(Bz+  denotes  benzylium,  Tr+  is  tropylium).  Channels  1 .3  and 
1.5— 1.7  are  minor  under  the  present  conditions  and  can  be 
mostly  ignored.  The  bulk  of  this  article  focuses  on  reactions 
1.1,  1.2,  and  1.4.  In  addition,  at  temperatures  approaching  600 
K,  collisionally  stabilized  C9Hi2+  ions  can  pyrolyze, 

C9H12+(+M)  —  C7H7+  +  C2H5(+M)  kg  (1.8) 

Secondary  reactions  occur  but  are  accounted  for  by  extrapolating 
to  zero  [C9H12]. 

There  is  also  a  practical  aspect  to  the  measurements.  The 
bottleneck  for  combustion  of  hydrocarbon  fuels  often  is  the 
ignition  delay  caused  by  the  primary  fragmentation  of  the  fuel 
molecules.  Injection  of  radicals  or  ions  can  considerably  speed 
up  the  ignition  and  may  be  required  for  advanced  hydrocarbon- 
fuelled  airbreathing  propulsion  systems.1  The  present  work 
contributes  to  this  practical  aspect  of  combustion  technology 
under  conditions  more  appropriate  to  real  life  engines. 

2.  Experimental  Technique 

Our  experiments  were  conducted  in  the  TIFT  which  has  been 
described  in  detail  before;4’16  therefore  the  technique  is  only 
briefly  summarized  here.  The  TIFT  is  similar  to  low-pressure 
flow  tubes  except  that  larger  flow  rates  are  used  to  achieve 
higher  pressures  and  Reynolds  numbers.  A  liquid  nitrogen 
storage  vessel  supplies  the  N2  carrier  gas,  which  is  preheated 
in  a  sidearm.  The  ions  are  created  upstream  in  an  off-axis  corona 
discharge  source  and  are  entrained  into  the  N2  flow  in  the 
sidearm.  The  gas  enters  the  flow  tube  where  the  neutral  reagent 
is  injected  through  a  moveable  on-axis  tube.  At  the  end  of  the 
flow  tube,  most  of  the  gas  is  removed  by  a  large  mechanical 
pump  while  a  small  fraction  is  sampled  through  a  150  jum  orifice 
in  a  truncated  nose  cone.  The  core  of  the  ensuing  supersonic 
expansion  is  sampled  through  a  skimmer  into  a  mass  spectrom¬ 
eter.  The  ions  are  analyzed  by  a  quadrupole  mass  filter  and 
counted  by  a  discrete  dynode  electron  multiplier. 

To  elevate  the  temperature,  six  zones  of  heating  are  used  and 
the  thermocouple  that  monitors  the  gas  temperature  is  main¬ 
tained  to  ±2  K.  The  main  flow  tube  is  heated  in  four  zones, 
two  short  zones  at  either  end,  a  long  middle  section,  and  a  zone 
inside  the  vacuum  chamber  just  upstream  of  the  nose  cone.  The 
connection  between  the  corona  discharge  tube  and  the  sidearm 


so  100 

Mass  (amu) 


Figure  1.  Typical  low-resolution  mass  spectrum  of  the  C>2+  +  C9H12 
system.  The  dashed  peak  at  mass  32  corresponds  to  the  undepleted 
02+  signal  before  the  addition  of  propylbenzene.  The  peak  at  mass 
133  corresponds  to  a  secondary  reaction  of  benzylium  (C?H7+)  with 
C9H12;  see  eq  2.1. 

is  also  heated.  Last,  it  was  found  that  preheating  of  the  buffer 
gas  is  needed.  During  the  course  of  this  work,  the  preheater 
design  was  changed  to  allow  for  higher  temperatures.  The  new 
preheater  consists  of  two  Mellen  split  type  (clam  shell)  tubular 
heating  elements  that  surround  a  tubular  piece  of  copper;  the 
gas  flows  through  narrow  channels  (which  traverse  the  copper 
tube  along  the  axis)  that  have  a  radial-spoke  cross  section.  The 
pattern  was  chosen  to  promote  efficient  heating. 

A  typical  mass  spectrum  is  shown  in  Figure  1,  with  and 
without  the  propylbenzene  reactant  gas.  In  the  absence  of 
propylbenzene,  the  02+  signal  is  dominant  (>99%). 

Common  impurities  are  H3O+(H2O)0,i  and  02+(H20).  Since 
the  relative  concentrations  of  the  impurities  were  usually  less 
than  2%  (sometimes  much  less),  no  corrections  to  the  recorded 
branching  ratios  were  made.  When  propylbenzene  is  added,  ions 
at  78,  91,  92,  105,  120,  and  133  amu  are  observed.  In  Figure  1, 
the  92  amu  peak  is  hidden  in  the  large  91  amu  peak.  Higher 
resolution  spectra  are  used  to  separate  the  two  peaks,  and  an 
isotope  correction  is  made.  The  mass  78  and  105  amu  peaks 
are  present  at  the  same  levels  as  impurities,  ca.  0.5—2%  at  low 
temperatures.  At  higher  temperatures,  their  yields  increase  and 
can  be  investigated.  The  dissociative  and  nondissociative  charge 
transfer  products  at  91  and  120  amu,  respectively,  dominate. 
The  peak  at  mass  133  is  due  to  a  reaction  between  benzylium 
and  propylbenzene3,  that  is, 

C7H7+(Bz+)  +  C9H12  -  C10H,3+  +  C6H6  (2.1) 

Branching  ratios  of  the  ion  products  were  recorded  as  a  function 
of  the  neutral  reagent  concentration.  Secondary  chemistry  was 
taken  into  account  by  extrapolating  to  zero  concentration.  The 
accessible  temperature  and  pressure  ranges  were  limited  by  the 
allowable  intensity  of  impurity  ions  (<2%)  at  low  T  and  high 
P,  thermal  decomposition  of  C9Hi2+  at  high  T,  and  small  signals 
at  low  P. 

Up  to  523  K,  the  pyrolysis  reaction  1 .8  could  be  neglected. 
In  this  case,  the  mass  120  versus  91  ratio,  which  we  call  S/D, 
can  be  symbolically  written  as 

S  _  [C9H12+]  _  *4[M] 

D  [C7H7+]  k2  +  k,  {22) 

Here  it  is  assumed  that  at  zero  pressure,  all  of  the  C9H]2+ 
dissociates,  but  in  section  3,  it  is  shown  that  the  S/D  values 
exhibit  a  nonzero  intercept  at  low  pressures,  indicating  that  some 
of  the  C9Hi2+  transfer  product  is  stable  at  zero  pressure.  Minor 
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channels  are  excluded  in  the  dissociation  signal.  At  T  >  573 
K,  S/D  values  vary  with  the  injector  position  and  therefore  time, 
indicating  that  thermal  decomposition  of  is  occurring. 

Our  first  attempt  to  determine  the  pyrolysis  rate  of  C9Hi2+ 
involved  direct  observation  of  the  depletion  of  C9Hi2+  as  a 
function  of  distance  by  using  NO+  or  CeH.6+  as  reagent  ions. 
Because  the  ionization  potentials  of  CgHg,  NO,  and  C9H12  are 
similar,17  this  avoided  the  unwanted  presence  of  highly  vibra- 
tionally  excited  C9H]2+.  Unfortunately,  this  approach  yielded 
uninterpretable  results  and  was  abandoned.  Explanations  for  the 
failure  include  interference  from  the  isomerization  of  the 
benzene  cation  to  the  five-membered-ring  species,  fulvene,18 
and  hydrogen  abstraction  reactions  when  NO+  was  used  as  the 
charge  transfer  agent. 

The  method  that  was  successful  involves  the  deconvolution 
of  the  effects  of  time  and  pressure.  An  advantage  over  the 
previous  method  is  that  (S/D)  in  the  absence  of  thermal 
decomposition  ( t  =  0),  called  (S/D)o,  can  also  be  obtained.  S/D 
was  measured  in  the  usual  way  at  several  distances  (i.e.,  reaction 
times)  at  the  same  T  and  P.  This  was  repeated  at  numerous 
pressures,  P  =  30—250  Torr,  at  573  and  603  K. 

A  steady-state  solution  of  the  kinetics,  modeled  by  eqs  1.1  — 
1.4  and  1.8  which  include  stabilization,  dissociation,  and  thermal 
dissociation,  yields 


*4[M] 

(k2  + 1 13) 


[1  -  exp(— V)] 


D 


V 1 1  + 


A4[M] 


(*2  + *3 )/  \(*2  +  *3). 


*4[M] 


(2.3) 


[1  -  exp(-y)] 


Minor  channels  are  neglected.  Extrapolating  to  zero  time  leads 
back  to  eq  2.2,  as  expected;  this  value  is  (S/D)o.  At  high 
temperatures  and  nonzero  time,  substitution  of  eq  2.2  into  eq 
2.3  yields  the  following: 


_ (S/D)0[l  -  exp(-y)] _ 

D  ~  V(  1  +  (S/D)0)  -  (S/D)0[l  -  exp(-V)] 


(2.4) 


We  assume  that  the  decrease  of  the  product  of  concentrations 
[C9Hi2][C>2+]  during  the  course  of  the  reaction  is  insignificant 
since  we  extrapolate  to  zero  extent  of  reaction.  Nonlinear  least- 
squares  fits  of  the  time  dependence  of  S/D  at  each  pressure  to 
eq  2.4  yield  both  k%  and  (S/D)o.  Note  that  this  involves  both  an 
extrapolation  to  [C9H12]  =  0  and  the  fit  to  the  time  dependence, 
thereby  introducing  noticeable  scatter  in  the  data.  The  fitted 
values  at  573  K  varied  by  a  factor  of  2.  The  scatter  at  603  K  is 
less  since  the  pyrolysis  rate  is  larger.  The  factor  of  2  scatter  is 
an  estimate  of  the  accuracy.  Figure  2  shows  a  typical  time 
dependence  of  S/D  at  T  —  573  and  603  K  and  at  a  nitrogen 
pressure  of  100  Torr.  The  pyrolysis  decay  is  easily  detected. 
We,  therefore,  have  the  interesting  situation  of  a  chemical 
activation  process  which  is  initiated  by  a  charge  transfer  reaction 
and  which,  after  collisional  stabilization  of  the  highly  excited 
products,  is  followed  by  thennal  dissociation  of  the  stabilized 
species.  The  intrinsic  kinetics  of  chemical  activation  and  thermal 
dissociation  is  connected  and  must  be  inteipreted  in  a  consistent 
way,  and  this  is  done  here. 

The  fragmentation  of  C9Hi2+*  is  known  to  lead  to  two 
isomers  of  C7H7+,  Bz~  and  Tr+.15  The  former  contains  the  six- 
membered  aromatic  ring  and  is  reactive  with  propylbenzene  to 
produce  C9H1 1+  (mass  1 19)  and  CioHb+  (mass  133)  via  reaction 
2.1.  By  contrast,  the  latter  is  comprised  of  an  unreactive  seven- 
membered  ring.  This  difference  in  reactivity  allows  for  the  Tr+ 


t/s 

Figure  2.  Thennal  decomposition  of  CsH]2+.  [C9Hi2+]/[C7H7+]  is 
plotted  vs  reaction  time  at  pressures  of  100  Torr  of  N2.  Best  fits  to  eq 
2.4,  shown  as  curves,  yield  values  of  (S/D)o,  the  ratio  of  stabilization 
to  dissociation  in  the  absence  of  thermal  decomposition,  and  k%,  the 
first-order  thermal  decomposition  rate  coefficient  of  C9H]2+. 


Figure  3.  Rate  constants  k\  for  the  reaction  of  02+  with  GjHn  as  a 
function  of  pressure  at  various  temperatures. 

yield  to  be  measured  by  alternately  using  low  and  high 
propylbenzene  flows.  Low  concentrations  of  propylbenzene  are 
used  to  measure  the  total  relative  concentration  of  C7H7+, 
whereas  at  high  concentrations  Bz+  is  completely  removed. 
Subtraction  of  the  two  limits  yields  the  amount  of  Bz+.  At  423 
K,  the  Tr+  yield  was  found  to  be  about  1  %,  and  at  523  K  the 
value  was  roughly  7%.  At  573  K,  a  yield  around  20%  was 
found.  The  uncertainty  is  considerable  since  the  propylbenzene 
cation  starts  to  decompose  at  this  temperature  and  the  kinetics 
of  reaction  2.1  has  not  been  measured  at  high  temperature. 

3.  Results 

The  experimental  results  are  of  three  types:  (1)  the  total  rate 
constant  k\  of  the  charge  transfer  process  1.1,  (2)  branching 
ratios  as  a  function  of  temperature  and  of  pressure  (mostly  the 
ratio  of  to  C7H7+  but  the  other  channels  are  also 

reported),  and  (3)  the  thermal  decomposition  rate  constants  at 
high  temperature  as  a  function  of  pressure.  Rate  constants  k\ 
are  plotted  in  Figure  3,  and  S/D  values  are  listed  in  Table  1. 
Table  2  lists  the  thennal  decomposition  rates,  and  Table  3  lists 
the  ratio  of  mass  91  to  the  minor  products  and  to  the  fraction 
of  mass  91  that  is  tropylium. 

The  rate  constants  shown  in  Figure  3  show  no  dependence 
of  k\  on  temperature  or  pressure  over  the  ranges  423—523  K 
and  30—200  Torr.  The  uncertainty  of  individual  rate  constants 
was  about  ±25%.  Part  of  the  scatter  is  because  most  of  the 
effort  was  focused  on  the  more  interesting  branching  information 
where  low  depletion  is  used.  Rate  constants  are  better  obtained 
under  high  depletion.  The  average  value  of  all  measurements 


V  *  * 

Stabilization  and  Dissociation  of  CgH^  Ions 


J.  Phys.  Chem.  A,  Vol.  108,  No.  45,  2004  9655 


TABLE  1:  Stabilization  vs  Dissociation  Yields  S/D  = 


[C9H12+]/[C7H7+]  in  N2‘ 


77K 

P/Torr 

m 

S/D 

Cl 

Cl 

-  423 

15 

0.31 

1.18 

0.59 

1.44 

30 

0.65 

1.61 

50 

1.15 

2.28 

75 

1.69 

3.18 

448 

15 

0.29 

1.07 

0.86 

0.516 

30 

0.62 

1.07 

50 

1.05 

1.43 

75 

1.58 

1.69 

100 

2.11 

1.94 

150 

3.19 

2.22 

473 

15 

0.28 

0.87 

0.79 

0.323 

30 

0.60 

0.96 

50 

0.99 

1.13 

75 

1.51 

1.29 

100 

2.03 

1.48 

150 

3.04 

1.74 

523 

15 

0.26 

0.75 

0.68 

0.134 

30 

0.54 

0.66 

50 

0.92 

0.83 

75 

1.37 

0.86 

100 

1.81 

1.00 

573 

30 

0.505 

0.63 

0.54 

0.147 

50 

0.835 

0.67 

75 

1.26 

0.82 

100 

1.68 

0.90 

150 

2.52 

0.95 

200 

3.35 

1.08 

250 

4.19 

1.36 

603 

30 

0.474 

0.46 

0.48 

0.0599 

50 

0.798 

0.54 

75 

1.19 

0.65 

100 

1.59 

0.72 

150 

2.39 

0.66 

200 

3.19 

0.92 

250 

4.00 

1.32 

0  [N2]  in  1018  molecule  cm  3, 

c2  in  10-18 

molecule-1 

cm3. 

TABLE  2: 

Thermal  Dissociation  Rate  Constants  ks  for 

c9h12+  - 

c7h7+ 

+  c2h5 

77K 

P/Torr 

fa/s~' 

77 K 

P/Torr 

kg/  s  1 

573 

30 

75 

603 

30 

101 

50 

37 

50 

109 

75 

44 

75 

109 

100 

32 

100 

98 

150 

22 

150 

75 

200 

27 

200 

109 

250 

34 

250 

138 

was  found  to  be  k\  =  2.1  x  10  9  cm3  molecule-1  s-1,  in 
agreement  with  our  previous  measurements,1,3  and  confirms  the 
accuracy  of  the  TIFT  approach.  This  value  is  also  in  good 
agreement  with  the  Langevin  rate  constant  for  capture  of  C9H12 
by  O2"1",  which  is  calculated  as  =  1.89  x  10-9  cm3  molecule-1 
s-1  on  the  basis  of  an  estimated19  polarizability  of  C9H12  of  a 
=  16.1  x  10-24  cm3  (contributions  from  polarizability  anisotro¬ 
pies  and  quadrupole  and  permanent  dipole  moments  at  this  level 
can  be  neglected). 

The  measurements  of  S/D  (or  (S/D)0  at  higher  temperatures) 
are  presented  in  Table  1  and  plotted  as  a  function  of  the 
temperature  and  the  buffer  gas  concentration  in  Figure  4.  The 
data  show  that  lower  temperatures  and  higher  pressures  favor 
the  formation  of  the  nondissociative  product,  that  is,  give  larger 
S/D.  The  measurements  exhibit  a  roughly  linear  dependence 
on  [N2]  or  equivalently  pressure.  One  also  notices  an  apparent 
intercept  of  the  S/D  versus  [N2]  curves  as  [N2]  — * -  0  which  is  in 
contrast  to  the  predictions  from  eq  2.2. 

The  present  results  are  similar  to  our  previous  observations 
for  the  02+  +  ethylbenzene  reaction.  The  intercept  indicates 
that  there  is  a  fraction,  a,  of  reaction  1.1  that  leads  to 


TABLE  3:  Ratio  of  Concentrations  of  Major-  to  Minor-Ion 
Fragments  from  the  Reaction  02+  +  C9Hi2  ~ *  C7H7+  + 
Other  (Minor)  Ion  Fragments 


77  K 

P/Torr 

[C7H,+]/ 

[C7H8+] 

[C7H7+]/ 

[C6H6+] 

[C7H7+]/ 

[CsH9+] 

[Tr+]/ 

[C7H7+]tot0,4 

423 

15 

17.8 

23.0 

51.0 

0.010 

30 

16.3 

25.9 

64.8 

50 

18.2 

28.4 

78.0 

75 

39.8 

60.0 

448 

15 

16.1 

30 

15.2 

20.3 

50 

16.7 

23.1 

83.0 

75 

17.3 

26.9 

75.4 

100 

19.6 

26.5 

69.0 

150 

18.3 

28.8 

79.3 

473 

15 

25.0 

40.4 

0.010 

30 

18.5 

22.0 

63.3 

50 

14.8 

24.8 

67.3 

75 

18.8 

25.8 

87.8 

0.041 

100 

21.3 

25.4 

67.8 

0.030 

150 

16.4 

93.5 

200 

16.1 

523 

15 

16.4 

19.1 

25.2 

0.032 

30 

16.3 

21.9 

53.6 

50 

17.8 

20.8 

60.1 

0.091 

75 

18.3 

25.5 

76.4 

100 

18.9 

23.7 

62.3 

150 

25.8 

20.0 

53.6 

0.070 

200 

21.1 

18.4 

65.9 

“  Fraction  of  Tr+,  the  seven-membered-ring  tropylium  isomer  of 
C7H7+;  the  aromatic  benzylium  isomer  is  the  major  species.  b  Data  of 
[Tr+]/[C7H7+]tot  at  573  K,  about  18%  and  invariant  from  30  to  200 
Torr,  are  not  included  here  due  to  interference  from  thermal  decomposi¬ 
tion. 


Figure  4.  Ratio  of  stabilization  to  decomposition  vs  [N2].  Temperatures 
range  from  423  to  603  K,  while  pressure  was  varied  from  about  15  to 
250  Torr.  At  T>  573  K,  thermal  decomposition  is  considerable  and  is 
taken  into  account. 

nondissociative  products.  Therefore,  the  mechanism  of  reactions 
1.1  — 1.8  is  extended  to  (excluding  minor  channels) 


o2+  +  c9h12-o2  +  c9h12+ 

O^cap 

(3.1) 

o2+  +  c9h12  —  o2  +  c9h12+* 

(1  -  a)*cap 

(3.2) 

c9h12+*  —  c7h7+  +  c2h5 

m 

(3.3) 

C9H12+*  +  M  —  C9HI2+  +  M 

(3.4) 

C9H12+  (+  M)  —  CvH7+  +  C2H5  (+  M)  itg 

(3.5) 

C7H7+  in  reaction  3.3  stands  for  the  sum  of  Bz+  and  Tr+  but, 
as  shown  in  Table  3,  is  mostly  the  former.  E  denotes  the  internal 
energy  of  the  chemically  activated  propylbenzene  cations,  and 
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TABLE  4:  Energy  Transfer  Properties  for  Highly  Excited 
Propylbenzene  Ions" 


77K 

k(E)/yc  s  1 

7(7) 

7c  (T=  0) 

— 1 JAE\/hc  cm  : 

423 

7.18  x  10s 

0.278 

0.195 

336 

448 

2.34  x  109 

0.240 

0.069 

119 

473 

1.12  x  109 

0.207 

0.168 

288 

523 

2.89  x  109 

0.152 

0.089 

153 

573 

2.87  x  109 

0.110 

0.123 

212 

603 

7.33  x  109 

0.0905 

0.059 

101 

a  See  the  text;  analysis  with  eqs  4.4  and  4.5;  E^Jhc  =  27  030  cm-1, 
k(Ect,)  —  3.90  x  107  s-1,  s*  =  7.608,  Z  =  6.5  x  1 0— 10  cm3  molecule"1 
s_l,  ci  and  c2  from  Table  1. 


k(E)  is  the  specific  rate  constant  for  their  fragmentation.  We 
identify  Z  with  the  Langevin  capture  rate  constant  between 
C9Hi2+  and  M:  Z  =  6.5  x  10-10  cm3  molecule-1  s-1.  S/D  is 
represented  in  this  mechanism  as 

S_[C9H12+]  q  ycz[N2] 

D  [C7H7+]^l-a  (1  ~  a)k(E) 

Linear  extrapolations  of  S/D  versus  [N2]  in  Figure  4  to  [N2]  — • - 
0  yield  intercepts  on  the  order  of  0.66  for  all  temperatures  giving 
a  =  0.40,  similar  to  the  value  of  a  =  0.44  found  in  our  02+  + 
CgHio  experiments.7 

The  apparent  intercepts  in  Figure  4  are  significantly  larger 
than  S/D  values  observed  in  our  previous  measurements1  of  the 
02+  +  CgHq  reaction  at  lower  pressures;  for  example,  measure¬ 
ments  in  the  selected  ion  flow  tube  and  the  high-temperature 
flowing  afterglow  at  0.5  Torr  of  N2  at  500  K  gave  S/D  =  0.15. 
Furthermore,  experiments  from  a  single-collision  guided-ion 
beam  apparatus  led  to  values  of  S/D  =  0.19  and  0.1 1  at  300 
and  650  K,  respectively.  For  CgHio4-,  we  attributed  this 
“nondissociative”  fraction  to  incomplete  energy  transfer  in  some 
of  the  charge  transfer  encounters,  which  is  quenched  at  lower 
pressures  than  studied  here.  Several  types  of  charge  transfer 
processes  were  postulated  and  were  analyzed  in  detail  in  ref  7. 
The  present  work  employs  an  analogous  interpretation.  The  three 
processes  are  (1)  resonant  charge  transfer  in  which  the  difference 
of  the  ionization  energies  of  02  and  C9H12  into  C9FIi2+  is 
transferred,  (2)  complex  formation  where  the  energy  transfer  is 
statistical,  and  (3)  resonant  charge  transfer  where  the  02  neutral 
is  produced  in  an  excited  state. 

The  results  from  Table  1  were  fit  to  a  linear  expression 
according  to  eq  3.5,  S/D  wq  +  c2[N2],  and  the  fit  parameters 
ci  and  c2  are  summarized  in  Table  1.  The  values  for  k(E)/yc  = 
Z(1  +  ci)/c2  are  given  in  Table  4.  k(E)/yc  in  this  representation 
corresponds  to  averages  over  suitable  energy  distributions  of 
C9H12"1"*.  A  more  detailed  analysis  is  given  below. 

.  The  rate  constant  for  the  thermal  dissociation  of  C9Hi2+  has 
been  measured  at  573  and  603  K.  The  average  values  are  fa  = 
33  s-1  at  573  K  and  fa  =  106  s-!  at  603  K  and  do  not  depend 
on  pressure,  indicative  of  the  high-pressure  limit.  Our  analysis 
presented  below  is  in  accord  with  this  conclusion.  The  corre¬ 
sponding  Arrhenius  expression  is  fa  ~  5.7  x  1011  exp(— 112 
kJ  mol-1//? 7)  s-1.  The  factor  of  2  uncertainty  in  the  rates  and 
the  small  temperature  range  lead  to  uncertainty  on  the  order  of 
70  kJ  mol-1  in  the  activation  energy.  If  one  assumes  the 
activation  energy  to  be  the  dissociation  energy,  E()  =  167  kJ 
mol-1,15  the  two  values  of  fa  correspond  to 

ks  «  2.1  (±0.9)  x  1016  exp(— 167  kJ  moF1//??)  s_1  (3.7) 

While  this  expression  appears  more  reasonable,  an  internally 
consistent  derivation  of  fa  is  given  in  the  discussion  section. 


All  minor  channels  decrease  with  increasing  pressure  (relative 
to  [C9Hi2+]).  However,  the  data  had  too  much  scatter  to  ascertain 
if  individual  channels  had  different  dependences.  No  attempt 
was  made  to  correct  these  channels  for  impurities,  since  the 
procedure  was  uncertain.  The  most  abundant  of  the  minor 
channels  is  C-jHfa ,  which  is  approximately  5—6%  of  the  C7H7+ 
channel.  There  is  no  temperature  dependence  within  our 
uncertainty,  as  is  true  for  all  minor  species.  CgHf,4  and  CgH9+ 
are  about  4  and  2%  of  the  C7H7+  signal,  respectively.  The  Tr+ 
is  small  except  in  the  temperature  range  where  thermal  dis¬ 
sociation  becomes  important,  where  20%  of  C7H7+  was 
observed  as  tropylium.  However,  the  monitor  ion  chemistry  was 
not  measured  at  high  temperature. 

4.  Analysis  of  Experimental  Results 

4.1.  Energy  Distributions  from  Charge  Transfer.  The 
observation  of  S/D  values  between  0. 1  and  0.2  under  single¬ 
collision  conditions  indicates  that  part  of  the  charge  transfer 
reaction  1.1  generates  C9Hi2+*  with  very  low  energy.  Either 
the  energy  of  this  fraction  is  below  the  dissociation  threshold 
or  can  be  stabilized  by  radiative  cooling  in  ~10-2  s.  We 
speculate  that  this  fraction  is  due  to  production  of  electronically 
excited  02.7  The  intercepts  in  Figure  4  are  larger  than  the  low- 
pressure  limiting  value,  providing  evidence  for  a  second  channel 
where  only  partial  energy  transfer  occurs.  We  suggest  that  this 
involves  fonnation  of  the  complex  (02-C9H]2)+  which  upon 
dissociation  distributes  energy  statistically  between  02  and 
C9Hi2+.  Finally,  near-resonant  charge  transfer  deposits  the 
difference  between  the  ionization  energies  of  02  and  C9Hj2  into 
vibrational  energy  of  C9H12+;  at  the  same  time,  the  thermal 
energies  of  02+  and  C9Hj2  are  assumed  to  be  nearly  retained 
as  internal  energies  of  02  and  C9Hi2+,  respectively.  Mainly  the 
resonant  channel  is  the  one  affected  by  pressure  in  the  present 
experiments.  In  summary,  the  three  channels  are 

02+  +  C9H12-02  +  C9H12+*  ga  (4. 1  a) 

02+  +  C9H12  -  (02-C9H12)+  -o2  +  C9H12+*  gb 

(4-lb) 

02+  +  C9H12  -  02*  ('A  or  ’2)  +  C9H12+  gc  (4.  lc) 

with  branching  fractions  ga  for  near-resonant  charge  transfer, 
gb  for  complex-forming  charge  transfer,  and  gc  for  near-resonant 
or  complex-fonning  charge  transfer  with  electronic  excitation 
of  02,  such  that  ga  +  gb  +  gc  =  1- 

We  modeled  the  energy  distributions  g(E,T)  of  C9Hj2+  arising 
from  channels  4.1a— c  by  analogy  to  ref  7.  For  convenience, 
the  thermal  vibrational  distribution  was  calculated  with  the 
density  functional  theory  (DFT)  frequencies  of  C9H]2+  derived 
in  ref  15.  Employing  the  energy  parameters  from  the  same 
reference,  the  difference  of  the  ionization  energies  of  02  and 
C9H]2  amounts  to  Ecfahc  =  27  030  cm-1  whereas  the  dissocia¬ 
tion  energy  of  C9H1F  into  C2H5  +  C7H7+  (benzylium)  at  0  K 
is  estimated  to  be  Eo/hc  =  13  950  cm-1. 

This  value  will  be  of  particular  importance  for  the  analysis 
of  fa.  Figure  5  shows  a  representative  modeling  result  for  g(E, 
473  K)  using  ga  =  0.60,  gb  =  0.27,  and  gc  =  0. 13  and  assuming 
that  the  (02— C9H]2)+  complex  distributes  its  energy  statistically 
over  the  vibrations  of  02  and  C9Hi2+  and  the  translation  and 
rotations  of  the  relative  motion. 

Unfortunately,  the  available  data  on  S/D  at  low  and  inter¬ 
mediate  pressure  are  not  sufficient  to  fix  the  branching  ratios 
ga  and  gb  with  certainty.  This  prevents  inclusion  of  a  temperature 
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Figure  5.  Energy  distributions  g(E,T)  (in  1/crrr1)  of  excited  C9Hi2+* 
from  the  charge  transfer  C>2+  +  n-CgHn  (T  —  473  K,  representative 
calculations  with  assumed  branching  ratios  ga  =  0.60  (curve  at  higher 
energies)  and  gt,  =  0.27  (curve  at  lower  energies),  see  the  text). 

dependence.  The  value  of  gc  ~  0.13  follows  from  the  single 
collision  results.  For  this  reason,  we  cannot  make  an  analysis 
which  accounts  for  the  full  energy  distribution  g(E,T).  Instead, 
we  have  to  limit  ourselves  to  an  analysis  based  on  eq  3.6.  It 
was  shown  in  our  C>2+  +  QHio  study  that  this  simplified 
analysis  leads  to  an  overall  uncertainty  of  about  a  factor  of  2  in 
%k£i.  The  fraction  1  —  a  of  channel  3.2  includes  not  only  ga 
from  channel  4.1a  but  also  the  high-energy  portion  from  channel 
4.1b.  Figure  5  illustrates  the  corresponding  overlap  of  the 
distributions  ga  and  gb.  Isomerization  was  considered  as  a  source 
of  the  intercepts  but  was  dismissed  because  the  temperature 
dependence  of  the  intercepts  would  be  much  larger  for  a  process 
with  a  large  barrier. 

4.2.  Specific  Rate  Constants  for  the  Dissociation  of 

C<)Hi2+*.  As  the  product  Z°/A£i  is  derived  from  a  relative  rate 
measurement,  its  precision  depends  on  that  of  the  reference  rate 
constant,  that  is,  the  specific  rate  constant  k(E)  for  dissociation 
of  Cc>Hi2+*,  in  addition  to  uncertainties  of  the  energy  distribution 
g(E,T).  Unfortunately,  there  exists  experimental  information  on 
HE)  only  from  a  single  study.15  These  data  have  been  fitted  by 
a  rigid  activated-complex  RRKM  model  in  ref  15,  the  results 
of  which  can  be  represented  by  the  expression 

HE)  =  3.6  x  108  s_1  [(£  -  E0)/hc  18310  cm-1]7'776  (4.2) 

with  EH  he  =  13  950  cm-1.  The  experimental  data  in  the  range 
(1  —4)  x  108  s"1,  of  course,  are  well  reproduced  by  eq  4.2  since 
they  defined  the  fitted  activated  complex  frequencies  of  the 
RRKM  representation. 

The  k(E)  data  from  ref  1 5  were  measured  at  energies  not  much 
above  those  reached  in  the  present  experiments.  They  are, 
therefore,  sufficient  to  evaluate  our  measured  S/D  values  with 
respect  to  energy  transfer.  However,  without  further  consider¬ 
ations  they  cannot  be  used  to  establish  a  relation  to  the  thermal 
dissociation  rate  constants.  First,  there  is  the  lack  of  a  proper 
account  for  the  /-dependence  of  the  specific  rate  constants, 
HE,J),  and  in  particular  of  the  threshold  energies  Eq(J),  which 
are  essential  for  the  calculation  of  the  thermal  rate  constant, 
H.  Second,  the  RRKM  fit  corresponds  to  a  vibrationally 
adiabatic  treatment  of  the  dynamics.  However,  it  was  shown  in 
ref  20  that  rigid  activated  complex  RRKM  theory  does  not 
provide  an  adequate  approach  to  k(E)  for  simple  bond  fission 
processes  when  the  Massey  parameter  for  the  transitional  mode 
dynamics  corresponds  to  nonadiabatic  behavior.  This  is  par¬ 
ticularly  true  for  systems  with  transition  state  shifting,  taking 
place  on  a  potential  energy  surface  with  both  long-range 
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E  /  he  cm'1 

Figure  6.  Specific  rate  constants  HE)  for  the  dissociation  C<)Hi2+  — • 
C7H7+  +  C2H5  (experimental  points  from  ref  15;  however,  see  ref 
15a):  long— short  dashed  curve,  PST;  dashed  curve,  SACM/CT 
modeling  from  ref  20;  dotted  curve,  RRKM  fit  from  ref  15;  full  line, 
empirical  representation  by  eq  4.3.  All  curves  (except  PST)  are  fine- 
tuned  to  the  experimental  HE)  at  E/hc  =  32  260  cm-1  from  ref  15;  see 
the  text. 

electrostatic  and  short-range  valence  character.  The  Massey 
parameter  for  the  CgH]2+  dissociation  indeed  signals  nonadia¬ 
batic  dynamics20  so  that  an  RRKM  treatment  is  good  only  for 
fitting  purposes.  There  is  an  interesting  apparent  contradiction: 
on  one  hand,  kinetic  energy  distributions  of  the  fragments  often 
are  well  represented  by  loose  activated  complex  transition  state, 
or  phase  space,  theory  (PST);  on  the  other  hand,  specific  rate 
constants  k(E)  are  far  below  PST  predictions.  This  apparent 
contradiction  can  be  resolved  by  realizing  the  nonadiabatic 
character  of  the  dynamics:  (1)  the  dissociation  bottleneck 
governing  k(E)  is  located  in  the  range  of  the  potential  where 
dynamical  constraints  arise  from  the  marked  anisotropy  of  the 
potential  such  that  k(E)  is  far  below  PST  values;  (2)  nonadiabatic 
dynamics  changes  the  kinetic  energy  distribution  among  the 
fragments  beyond  the  bottleneck  far  into  the  range  where  the 
potential  becomes  increasingly  isotropic  such  that  PST  applies. 

In  the  modeling  study  of  ref  20,  we  have  designed  long-range/ 
short-range  switching  potentials  and  calculated  k{EJ)  using  the 
statistical  adiabatic  channel  model/classical  trajectory  (SACM/ 
CT)  approach.  Our  results  for  ethylbenzene  cation  dissociation 
were  found  to  be  in  fairly  good  agreement  with  experimental 
results  of  k(E)  over  the  range  103— 108  s~'.  This  provides 
confidence  in  the  k{E)  for  the  dissociation  of  n-propylbenzene 
cations  which  was  also  modeled  in  ref  20.  These  modeling 
results  could  well  be  represented  by  the  empirical  expression 

k(E,J=  0)  = 

3.6  x  108  s_I  [(£  -  E0(J—0))/hc  18310  cm-1]6'541  (4.3) 

Equation  4.3  is  fine-tuned  to  the  experimental  values  k(E)  = 
3.6  (±0.4)  x  108  s_1  at  E/hc  =  32  260  cm"1.  One  finds  that 
HE)  is  somewhat  closer  to  PST  results  at  lower  energies  and  it 
falls  increasingly  below  PST  at  higher  energies.  This  behavior 
is  illustrated  in  Figure  6  where  k(EH=  0)  from  PST  is  compared 
with  the  RRKM  fit  shown  in  eq  4.2,  the  empirical  expression 
of  eq  4.3,  the  SACM/CT  results,20  and  the  experimental  points.15 
Equation  4.3  has  a  similar  energy  dependence  to  that  for  CgHjo+ 
dissociation  where  an  exponent  of  4.824  was  found.  The  larger 
number  of  medium-  to  low-frequency  oscillators  in  C9H12"  is 
responsible  for  the  larger  exponent  in  eq  4.3. 

CgHio+  and  CgHi2+  have  nearly  equal  dissociation  energies. 
Because  of  the  larger  number  of  oscillators,  one  would  expect 
that  k(E),  at  the  same  E  —  Eq,  is  smaller  for  C9Hi2+  than  for 
C8Hio+.  However,  the  opposite  is  true:  for  the  range  of  interest 
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in  the  present  work,  k(E)  apparently  is  about  a  factor  of  2  larger 
for  C9Hi2+  than  the  corresponding  value  for  C8Hio+.  Again, 
the  reason  is  that  there  are  more  medium-  to  low-frequency 
oscillators  in  C9Hi2+.  As  a  consequence,  higher  pressures  for 
collisional  stabilization  are  required  in  C9Hi2+  than  in  CgHio+. 
This  is  illustrated  by  comparing  Figure  4  of  the  present  work 
with  Figure  1  of  ref  7. 

4.3.  Collisional  Energy  Parameters  for  Excited  C()Hn+*. 
An  evaluation  of  S/D  versus  [M]  was  made  in  two  stages  as 
we  did  previously.7  First  it  was  assumed  that  the  nonresonant 
charge  transfer  channels  4.1b  and  4.1c  only  contribute  to  the 
apparent  intercepts,  and  second,  the  influence  of  the  energy 
distribution  of  channel  4.1b  was  taken  into  account.  The  former, 
a  simplified  analysis,  was  preferred  in  the  present  work. 
Nevertheless,  we  also  show  the  results  from  a  representative 
simulation  of  the  more  precise  analysis.  However,  we  are  unable 
to  fix  the  required  parameters  ga  and  gb  with  the  same  certainty 
that  was  possible  for  CsHio+.  The  analysis7  from  eq  3.6  leads 
to  the  expression 

sk^  +  _3mkz^_ 

D  1-0  (1  -  a)«£t, 


Figure  7.  Modeling  of  the  dependence  of  S/D  on  [N2]  (energy 
distributions  g(E,T)  from  section  4.1  with  representative  assumed  values 
ga  =  0.60,  gb  =  0.27,  and  gc  =  0.13;  experimental  points  as  in 
Figure  4). 

using  eqs  4.4— 4.6  and  accept  an  uncertainty  of  the  derived 
values  of  °AEi  of  about  a  factor  of  2. 

The  mean  value  of  °/AE\  from  Table  4  is 


with 


°/AE\/hc  =  -200  (±  1 00)  cm”1  (4.7) 


I(T)  = 


I  Ech-E0  yiEth  +  /_5j,\ 

ta-Uo  +  sJ  \  Ech  )  Qvib  eXp\  klj 


dEth 

(4.5) 


In  this  expression,  we  assume  HE)  to  be  of  the  form  k(E)  X  (E 
—  Eo)s*-1  used  in  eq  4.3,  p(El\l)  is  the  vibrational  density  of 
states,  and  2Vib  is  the  vibrational  partition  function  of  C9Hi2+. 
We  also  assume  that  °/AEi  is  proportional  to  E  and  does  not 
depend  on  T.  The  collision  efficiency  yc  according  to  ref  10 
can  be  related  to  °/AEl  by  the  expression 


Y c  ^  -°/AE\s* 
1  -  ycz  ~  E  -  E0 


(4.6) 


If  yc  1,  as  found  here,  the  left-hand  side  approaches  yc.  Table 
4  includes  the  corresponding  values  for  yc  (0  K)  ~  — °/AEls*/ 
(Ech  —  Eo)  and  I(T).  The  values  for  °/AEl ,  such  as  derived  from 
the  experimental  values  for  k(E)/yc,  are  also  given  in  this  table. 

The  detailed  analysis,  accounting  for  near-resonant  (eq  4.1a) 
and  nonresonant  (eqs  4.1b,c)  charge  transfer,  is  more  involved. 
It  has  been  described  in  ref  7  and  will  not  be  explained  here 
again.  We  only  illustrate  the  dependence  of  S/D  on  [M]  by 
employing  g(E,7)  from  section  4. 1  with  representative  branching 
ratios  ofga «  0.60,  gb  ~  0.27,  andgc  ~  0.13.  Figure  7  compares 
this  modeling  with  the  experimental  S/D  versus  [N2]. 

One  realizes  that  there  are  some  deviations  of  the  experi¬ 
mental  results  from  the  modeled  curves.  We  were  unable  to 
trace  the  origin  of  these  deviations,  in  particular  for  the 
measurements  at  423  K  which  are  not  included  in  Figure  7.  In 
addition,  the  modeled  S/D  versus  [N2]  plots  are  more  curved  at 
low  [N2]  than  the  experimental  data.  The  latter  may  indicate 
inadequate  branching  ratios  ga  and  gb  or  problems  with  the 
statistical  modeling  of  the  charge  transfer  channel  4.1b.  We 
varied  ga  and  gb  within  this  model  but  could  not  obtain  much 
better  agreement  with  the  experiments.  (Increasing  °/AEl  or 
decreasing  ga  would  have  produced  steeper  curves  of  S/D  versus 
[N2].)  Therefore,  we  limited  ourselves  to  the  simplified  analysis 


at  an  average  excitation  energy  of  E/hc  &  (30  000—33  000) 
cm-1.  This  value  of  °/AEi  is  in  good  agreement  with  the 
corresponding  result  for  stabilization  of  CrHio+*  in  collisions 
with  N2  which,  at  a  similar  excitation  energy,  gave  °/AE\/hc  = 
—285  (±150)  cm-1.  Although  there  remains  an  uncertainty  of 
about  a  factor  of  2  in  °/AEi,  we  find  again  that  %kEi  values  for 
excited  ions  are  similar  to  the  corresponding  °/AEi  values  for 
excited  neutral  molecules.21-25  In  contrast,  the  collision  frequen¬ 
cies  Z  for  ions  and  neutrals  are  markedly  different. 

4.4.  Relation  between  Specific  Rate  Constants  k(E)  and 
Thermal  Rate  Constants  k(T)  for  C9ll12+  Dissociation. 
Determining  both  specific  rate  constants  k(E)  and  thermal  rate 
constants  k%  for  the  dissociation  C9Hi2+  — *•  C7H7+(Bz)  +  C2H5 
offers  the  rare  possibility  to  check  the  internal  consistency  of 
the  derived  data.  This  analysis  requires  the  full  details  of 
unimolecular  rate  theory.  In  particular,  the  contributions  of 
angular  momentum  (quantum  number  J)  in  k(E,J)  need  to  be 
accounted  for,  both  in  energy-specific  and  in  thermal  measure¬ 
ments.  This  analysis  has  been  done  in  our  parallel  SACM/CT 
study20  where  k(E,J)  was  derived  from  trajectory  calculations 
employing  a  model  potential  which  was  fine-tuned  to  reproduce 
the  experimental  measurements  of  HE)-15  Based  on  these 
calculations  of  k(EJ),  thermally  averaged  rate  constants 
£rec,G(T)  for  the  recombination  C7H7+  +  C2H5  — *  C9Hi2+  were 
derived.  An  important  result  from  this  analysis  was  the 
demonstration  that  the  marked  difference  of  the  calculated  k(EJ) 
from  the  results  given  by  PST  also  occurs  in  the  thermal  rate 
constants:  k(EJ)  falls  increasingly  below  the  PST  results  with 
increasing  E,  and  krccfi(T)  analogously  falls  below  the  phase 
space  theory  results  with  increasing  T.  In  other  words,  k,QC  fi(7) 
is  expected  to  be  much  smaller  than  the  Langevin  rate  constant 
for  capture  of  C2H5  by  C7H7+.  On  the  basis  of  the  experimental 
k(E)  values,  the  recombination  rate  was  predicted20  to  be 

krcc{l(T)  ~  8.0  x  10”12(77600  K)”116  cm3  molecule”1  s”1 

(4.8) 

over  the  temperature  range  400—800  K.  In  contrast,  the 
Langevin  rate  constant  is  given  by  ki  =  1.0  x  10"9  cm3 
molecule-1  s-1. 


Stabilization  and  Dissociation  of  C9Hi2+  Ions 

To  compare  the  predicted  kKCy(T)  with  the  measured  fa  from 
section  3,  the  equilibrium  constant  was  calculated  as 

=  ([C7H7+]  [C2H5]/[C9H ,  2+])eq  =  kdJkrec  (4.9) 

and  k&iSSfi(T)  was  obtained  by  combining  Keq  with  kKCfi.  It  then 
has  to  be  decided  whether  the  measured  fa  is  less  than  the  high- 
pressure  limiting  dissociation  rate  constant  fobs, a  because  of 
falloff  effects.  For  internal  consistency,  we  calculated  Keq  with 
the  DFT  frequencies  from  ref  15,  keeping  in  mind  that  the 
largest  uncertainty  stems  from  the  bond  energy  Eo  which,  as 
before,  was  taken  to  be  EJhc  =  13  950  cm-1.  In  this  way, 
between  500  and  600  K,  we  obtained 


1.67  x  1027( 77600  K)~0'33  exp {-EJkT)  molecule  cm-3 

(4.10) 

To  assess  falloff  effects,  we  calculated  limiting  low-pressure 
rate  constants  fojSS.o  following  the  conventional  method  described 
in  ref  26.  The  relevant  quantities  entering  fo;SSio  are  Z  =  6.5  x 
10-10  cm3  molecule-1  s-1,  jic  &  0.2,  the  vibrational  density  of 
states  jOvib(£o)  =  5.4  x  1014/cm-1,  the  vibrational  partition 
function  Qvib(603  K)  =  8.2  x  105,  FE(600  K)  «  1.6,  and  Frol « 
F rot,max(603  K)  «  5.3.  Combining  /rdiss,ri  =  kKcfiKeq  from  eqs 
4.8—4.10  with  the  derived  fojSSio  indicates  that  the  center  of  the 
falloff  curve  foiSS,o([M])  =  foiss.fj  at  600  K  is  located  near 
[N2]center  «  4  x  1013  molecule  cm-3.  Falloff  curves  connecting 
foiss.o  with  /fdiss.a  have  been  illustrated  in  refs  26  and  27  for 
systems  of  comparable  size.  On  this  basis,  one  can  safely  assume 
that  the  measured  fo  corresponds  to  the  high-pressure  rate 
constant  foiSSia  within  a  few  percent.  The  absence  of  an 
observable  pressure  dependence  in  Table  2  confirms  this 
conclusion. 

Combining  the  modeled  faccfi  with  Keq  from  eqs  4.8-4.10 
predicts 

ks  =  1.34  x  10,6(77600  K)1'49  exp(-£0//c7)  s-1  (4.1 1) 

The  prediction  is  remarkably  close  to  our  direct  measurement 
of  fa  given  by  eq  3.7  with  E0  =  166.9  kJ  mol-1.  The  difference 
between  kKC&  and  foangevin  in  the  modeling  has  been  calibrated 
by  the  difference  between  the  measured  k(E)  and  k(E)  from  PST 
such  that  realistic  values  of  fo  should  have  been  produced.  The 
good  agreement  between  the  measured  fo  from  eq  3.7  and  the 
predicted  fo  from  eq  4.11  confirms  the  validity  and  internal 
consistency  of  the  described  approach.  The  agreement  between 
the  measured  and  the  modeled  value  also  supports  the  Eo  value 
of  Eo/hc  =  13  950  cm-1  (or  Eo  =  166.9  kJ  mol-1)  used  in  the 
analysis.  The  difference  of  a  factor  of  1 .6  in  fo  under  our  condi¬ 
tions  would  correspond  to  a  difference  of  about  190  cm-1  in 
Eo/hc.  The  energetics  of  the  dissociation  C9Hi2+  — ■  C7H7+  + 
C2H5,  therefore,  should  be  accurate  to  better  than  about  2.2  kJ 
mol-1. 

5.  Conclusions 

This  study  described  charge  transfer  generating  excited 
C9Hi2+  ions  which  then  dissociate  mainly  to  C7H7+  ions  and 
C2H5  radicals.  Calibrating  S/D  ratios  for  chemical  activation 
using  literature  values  of  the  specific  rate  constants  for  dis¬ 
sociation,  absolute  values  for  collisional  stabilization  rates  were 
obtained.  From  these  quantities,  average  energies  transferred 
per  collision  °/i!sE\/hc  =  —200  (±100)  cm-1  were  derived  for  a 
bath  gas  N2  at  excitation  energies  near  E/hc  =  32  000  cm-1. 
These  values  are  similar  to  those  measured  for  excited  CgHio+ 
ions,  where  °/AE\/hc  =  —285  (±150)  cm-1  was  obtained.  They 
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are  also  of  similar  magnitude  as  found  for  the  corresponding 
excited  neutral  molecules. 

The  conversion  of  specific  rate  constants  k(E)  into  high- 
pressure  thermal  dissociation  rate  constants  fojs s>r(7)  was  per¬ 
formed  by  unimolecular  rate  theory  in  the  form  of  SACM/CT 
calculations.  As  foiss,a(7)  has  a  particularly  large  temperature 
dependence,  these  quantities  also  allowed  us  to  confirm  the 
dissociation  energy  of  C9Hi2+  as  166.9  (±2.2)  kJ  mol-1  (at  0 
K).  The  given  analysis  provided  an  intrinsically  consistent 
description  of  rate  parameters  such  that  it  may  serve  as  an 
example  for  subsequent  studies  of  other  chemical/thermal 
activation  reactions  of  excited  molecular  ions. 
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